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It is shown that it is possible to quantitatively explain quantum Monte Carlo results for the 
Green's function of the two-dimensional Hubbard model in the weak to intermediate coupling regime. 
The analytic approach includes vertex corrections in a paramagnon-like self-energy. All parameters 
are determined self-consistently. This approach clearly shows that in two dimensions Fermi-liquid 
quasiparticles disappear in the paramagnetic state when the antiferromagnetic correlation length 
becomes larger than the electronic thermal de Broglie wavelength. 
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^] \ For more than thirty years, the Hubbard model has served as a theoretical basis for a wide range of interacting 
electron phenomena, from itinerant electron magnetism to, more recently, high-temperature superconductivity. In the 
latter context, one of the key issues has been the possibility of non- Fermi liquid behavior in the single-particle proper- 
ties. These properties have been intensively studied in the weak to intermediate coupling regime using paramagnon- 
type theories. |Q Q Such approaches are physically attractive since they describe the effect of low-lying collective 
modes on single-particle properties in a manner similar to that of phonons. However in these theories one faces the 
y—i \ difficulty that, contrary to the case of phonons, vertex corrections are not small (no Migdal theorem) and have to be 
£> ■ taken into account. 

In this letter, we introduce a simple approach that has the straightforward physical interpretation of paramagnon 
theories and uses no adjustable parameter. It does include vertex corrections, and achieves, for both two-particle and 
single-particle properties, better quantitative agreement with Monte Carlo simulations than previous theories || Q] 
|)j . Our analytical approach allows us to unambiguously interpret the Monte Carlo data and to show that in two 
I/"") | dimensions spin fluctuations destroy the Fermi liquid in the paramagnetic state when the antiferromagnetic correlation 
ON . length becomes larger than the thermal de Broglie wavelength of electrons. This corresponds to the disappearance 
of the quasiparticle peak and the appearance of a pseudogap, an issue that has been controversial mostly due to 
finite-size effects in Monte Carlo || . 

We consider the one-band Hubbard model on the square lattice with unit lattice spacing, on-site repulsion U and 
i-^j ■ nearest- neighbor hopping t. We work in units where the lattice spacing is unity, ks = 1 , h = 1 and t = 1. The theory 
has a simple structure that we now explain physically, postponing the formal derivation momentarily. 

The calculation proceeds in two steps: we first obtain spin and charge susceptibilities, then we inject them in 
the self-energy calculation. In the calculation of susceptibilities we make the approximation that spin and charge 
susceptibilities Xsp, Xch are given by RPA-like forms but with two different effective interactions U sp and U c h which 
are then determined self-consistently. The necessity to have two different effective interactions for spin and for charge 
is dictated by the Pauli exclusion principle = (n CT ) which implies that both Xsp and Xch are related to only one 
local pair correlation function (n^n^). Indeed, using the Fluctuation-Dissipation theorem in Matsubara formalism 
and the Pauli principle one can write: 

— V Xch »(?) = — V = n + 2 (-if (ntni ) - (1 - £) n 2 (1) 

^Vi + ^tWxo(g) 

where £ = for charge (c/i), I — 1 for spin (sp), (3 = 1/T, n = (rtf) + (rij), q = (q_,iq n ) with q the wave vectors of 
an N site lattice, iq n Matsubara frequencies and Xo(o) the susceptibility for non-interacting electrons. The value of 
(rif7ij_) may be obtained self-consistently fij by adding to the above set of equations the relation U sp = 3fj,(0) U with 

9n(°) = ( n 1 n l)/( n l)( n l)- 
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As shown in Rcf. M, the above procedure reproduces both Kanamori-Brueckner screening as well as the effect of 
Mcrmin- Wagner thermal fluctuations, giving a phase transition only at zero-temperature in two-dimensions. There is 
however a crossover temperature Tx below which the magnetic correlation length £ grows exponentially. Quantitative 
agreement with Monte Carlo simulations is obtained (tJ for all fillings and temperatures in the weak to intermediate 
coupling regime U < 8. 

We now turn to the discussion of the single-particle properties. In order, to be consistent with the two-particle 
correlation functions, the self-energy E CT (k) must satisfy the sum rule 



1 



A- «^E S " (*) G - e " fe " r = ^ ("T«l) . (2) 
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which follows from the definition of E CT (k) ||. Here, we encounter the same key quantity (n^ni) that appears in the 
sum rule for the susceptibilities Eq. (Q) . We find, using the formal approach described below, the following expression 
for E CT (jfe) 

E CT (fc) = Un_ a + j§ E PspXsvid) + U chX ch{q)\ G°(fc + q), (3) 
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which satisfies Eq. (H) with G<j replaced by G" on the left-hand side. This self-energy expression (|J) is physically 
appealing since, as expected from general skeleton diagrams, one of the vertices is the bare one U, while the other 
vertex is dressed and given by U sp or U c h depending on the type of fluctuation being exchanged. Eq.(|^) already 
gives good agreement with Monte Carlo data but the accuracy can be improved even further by requiring that the 
consistency condition (||) be satisfied with G CT instead of G° . To do so we replace U sp and U c h on the right-hand side 
of (^) by aU S p and aU c h with a determined self-consistently by Eq.(||). For U < 4, we have a < 1.15. This concludes 
the description of the structure of our theory. We now turn to the formal derivation. 

Formal aspects: We use the notation 1 = (Ri, ti), and an overbar (l) to indicate a space-time or spin variable which 
is integrated or summed over. Using the Matsubara Green's function G CT (1,3) in the presence of an external potential 
<f) a (2), one can write the exact equations for the generalized susceptibilities x , (1,3; 2) = 5G a (1,3) /8<j) a i (2) and 
the self-energy Y> a (1, 2) as 0: 

X ,,a- (1, 3; 2) = G a (1, 2) G CT (2, 3) 6^ + G CT (l, 6) G CT (7, 3) r£ ff (6, 7; 4, 5) (1, 5; 2) (4) 



E a (1, 2) = (1 - 2) - U Ga (l, 6) 1% (6, 2; 4, 5) (4, 5; 1) . (5) 

where T l £ a , = S~S a /SG a ' are irreducible vertices. Our approximation for spin and charge susceptibilities and for self- 
energy Eq. (||) is obtained from these exact equations by using as initial guess on the right-hand side of both equations 
a self-energy E^ and irreducible vertices which are functional derivatives of 

$ [G] = ^(U ch - U sp )n w (I) n w (T) + ^(U ch + U sp )n« (T) (l) (6) 

where n a (1) = G a (1, 1 + ). The initial self-energy E^ = 8^/8G a \4>=o is a constant and hence can be absorbed in 
a chemical potential shift so that all Green's functions on the right-hand side of Eqs. (Q) and (||) are bare ones 
G° (1,2). The corresponding particle-hole irreducible vertices r^°^ r = <5 2 $/((5G 2 )|0 = o are given in Fourier space by 
two constants 

Y^% = {U ch + U sp )/2 ■ T^ = (U ch -U sp )/2 (7) 

which are chosen so that the susceptibilities are consistent with the correct double-occupancy (n^n^) and with the 
exclusion principle (n 2 ) = (n a ), as described above. 

Important advantages of our iterative procedure are that: a) The approximation is conserving since E^ and 
are functional derivatives of a single functional <& |{J . b) The following identity 

E ff (l,3)G CT (3,2)=-t/(T T [n_ (T (l)^(l)V+(2)]) (8) 

is satisfied with E^ and G^ and it can be used as a convergence check, namely E^ (l, 3) (£> (3, 1+) « 

E^ (1,3) G [ f } (3,1+). For U < 4, this convergence criterion is satisfied at worse within 15%. The sum rule (0) 
follows from the identity (^) when 2 = 1 + . 

Comparisons with other theories and with Quantum Monte Carlo data: Fig. 1(a) shows G (k,r) for filling n = 0.875, 
temperature T = 0.25 and U = 4 for the wave vector on the 8x8 lattice which is closest to the Fermi surface, 
namely (tt, 0). For these parameters, size effects are negligible. Our theory is in agreement with Monte Carlo data 
and with the parquet approach but in this regime second-order perturbation theory for the self-energy gives the same 
result. This surprising performance of perturbation theory (see also pj) is a consequence of compensation between 
the renormalized vertices and susceptibilities (U sp < U, Xsp(q) > Xo(<z); U c h > U, Xchil) < Xo(o) )■ 
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Half-filling n — 1 is an ideal situation for numerical studies of low-energy phenomena since some of the allowed wave 
vectors on finite lattices lie exactly on the Fermi surface. Fig. 1(b) shows G(kp, t) for kp = (tt, 0) in a regime where the 
antiferromagnetic correlation length is growing exponentially. The Green's function is strongly renormalized since the 
non-interacting value is —0.5, independent of r. Our theory shows better agreement with Monte-Carlo than previous 
approaches. Despite the very large correlation length, finite-size effects for the Matsubara Green's function G(kp,-r) 
are not large and almost completely disappear when the system size becomes larger than the thermal de Broglie 
wavelength of the electrons (£th = ( v f) I^T) as we explain below. 

Our most dramatic numerical results addressing the issue of the Fermi liquid are shown in Fig. 2 where we plot 

z(T) = -2G(k F ,/3/2)= / p A fe"> ■ (9) 
v ' J 2?r cosh(/3o;/2) v ' 

This quantity z (T) is an average of the single-particle spectral weight A (kp, u>) within T = 1/(3 around the Fermi level 
(lu = 0) and it is a generalization of the usual zero-temperature quasiparticle renormalization factor z = 1/(1— dH/du)). 
For non-interacting particles z (T) is unity. For a normal Fermi liquid it becomes equal to a constant less than unity 
as the temperature decreases since the width of the quasiparticle peak scales as T 2 and hence lies within T of the 
Fermi level. The quantity z (T) is the best estimate of z one can get from Monte Carlo data for G(k, t). Moreover 
z (T) gives an estimate of A (kp, ui) around the Fermi surface even when the Fermi liquid does not exist. 

One can clearly see from Fig. 2 that while second-order perturbation theory exhibits typical Fermi-liquid behavior 
for z(T), both Monte Carlo data and a numerical evaluation of our expression for the self-energy lead to a rapid 
fall-off of z (T) below T x (for U = 4, Tx « 0.2 0). The decrease of z (T) in this regime is approximately z ~ T 2 
which clearly suggests non-Fermi liquid behavior. The physical origin of this effect is that the quasi-particles of the 
two-dimensional paramagnetic state become overdamped when the energy scale associated with the proximity to 
antifcrromagnetism SU = U m f,c — U sp {U m f,c = (Q>0)) becomes exponentially small or, more precisely, when 
the two-particle antiferromagnetic correlation length £ becomes larger than the single-particle thermal de Broglie 
wavelength = v f /ttT. 

Pseudogap: While size effects and statistical errors make continuation of the Monte Carlo data to real frequencies 
particularly difficult, in our approach we can make this continuation analytically to show that the above effect 
corresponds to a pseudogap. We first explain the behavior z ~ T 2 then do real frequency analysis. 

Since the spin susceptibility \sp 0) below Tx is almost singular at the antiferromagnetic wave vector Q = 
(tt, tt), the main contribution to E in Eq. (^) comes from iq n — and wave vectors (q — Q) 2 < near Q 

Approximating \s P (q, 0) in Eqs. ([!]) and @ by its asymptotic form \s P (q, 0)w2 U sp (, 2 (£~ 2 + (q — Q) 2 ; 
£q = 2xo(Q) 3 dl?^ anc ^ ^ = £° (U sp / 5U) x l 2 , the integrals over q can be done to obtain the asymptotic results 
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where 



l- 2 ^) no. 



e Ml) = u - - i UT m fc " + ^gZgg! + n. (ii) 

Here 1Z is the regular part which remains finite as T — * and a 2 = n — 2(n^n\) — C < 1 is the right- 
hand side of ([!]) minus corrections C that come from the sum over non-zero Matsubara frequencies (quantum 
effects) and from (q — Q) ^> . The corresponding value of z(T) can be written as the alternating series 
— 2G (kp, (3/2) — — 4T^^1 1 (—1)™/ (ik n + (J, — E (k F , ik n )). We can understand qualitatively both the tempera- 
ture and size dependence of the Monte Carlo data in Fig. 2 by using the first term of this series along with Eqs. ((TT 
and (ITT 



On the infinite lattice, £ starts growing exponentially below Tx, quickly becoming much larger than £ t h- This implies 
z(T) ~ T 2 . On finite lattices Eq.(0) gives £ ~ \/N below Tx which explains the size effect observed in Monte Carlo 
i.e. smaller z for smaller size N, (^th{Tx) ~ 5 for Fig. 2). 
The analytic continuation of E (kp, ik n ) in Eq. (|ll|) is 
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E«(k F ,c) = ^ 
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U + yfuj 2 + V 2 F £- 2 



- v/^ 2 + v 2 F i- 2 
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|2 + v F £- 2 
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(13) 



At small frequencies the regular part 1Z can be neglected so that exactly at the Fermi surface (uo = 0) the imaginary 
part of the self-energy for £ > £ t h increases exponentially when the temperature decreases, S"(kjr,o;) ~ f £/(£t/i£o) ~ 
exp (7T(T 2 £Q[/ S p/r) [Ly. This corresponds to a pseudogap: Instead of a quasiparticle peak, the spectral weight has a 
minimum at the Fermi level and two symmetrically located maxima away from it. Clearly this is not a Fermi liquid, 
yet the symmetry of the system remains unbroken at any finite T. By contrast jl2j with 3D, in 2D this effect exists in 
a wide temperature range T < Tx- These conclusions persist slightly away from half- filling. In particular, we do not 
find a quasiparticle peak in the pseudogap close to half-filling when £ ^> £th- This is different from the results inferred 
from a phenomenological zero- temperature calculation || (£ t h = oo) which physically corresponds to 1 -C £ <C ^th- 
in photoemission experiments on real quasi two-dimensional materials we predict a rapid decrease of the spectral 
weight and density of states at the Fermi level in a wide temperature range, from Tx to the Neel temperature T/v 
(T X -T N ~10 2 K), 
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Fig. 1. Comparison of our results for G(k, r) with Monte Carlo data, FLEX, parquet, and second-order pertubation 
theory, all on 8 x 8 mesh. Monte Carlo data and results for FLEX and parquet are from Ref. ||. a) n = 0.875, 
T = 0.25; b) n = 1, T = 0.17. 

Fig. 2. Temperature dependence of the generalized renormalization factor z defined in Eq.(|^). Lines are results of 
our calculations for infinite lattice and 16 x 16 mesh. Symbols are Monte Carlo data from Ref. H. 
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